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Abstract 

In the previous work we introduced a new static cylindrically symmetric vacuum solutions in 
Weyl coordinates in the context of the metric f (R) theories of gravity l|] . Now we obtain a 2- 
parameter family of exact solutions which contains cosmological constant and a new parameter 
as /3. This solution corresponds to a constant Ricci scalar. We proved that in f(R) gravity , 
the constant curvature solution in cylindrically symmetric cases is only one member of the most 
generalized Tian family in GR. We show that our constant curvature exact solution is applicable 
to the exterior of a string. Sensibility of stability under initial conditions is discussed. 
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I. INTRODUCTION 



Any phenomenal problem in theoretical physics commitment using of any thinkable 
method to attack. Acceleration of universe in now era is the most challenge for all branch of 
physics. If we are interested only on attribute it to a matter field, we envisage with different 
kind of models which we called them physical models or dark fluids. We present one of them 
as a simple phase transition in a known metric in [9]. Another approach is encroach to 
geometrical parts of Einstein field equations and working with a most general action, all 
these different proposals have been put forward with the hope of obtaining a consistent the- 
oretical background to the recent observation of an expanding universe which is seemingly 
not accessible through the Standard Model of cosmology. Modified or alternative theories 
of gravity is the paradigm under which all those theories which differ from the Einsteinian 
gravity are studied. One of these non-Einsteinian theories, that arose a lot of enthusiasm 
recently, is the so-called f (R) gravity in which a function, f (R), replaces the EinsteinHilbert 
(gravitational) Lagrangian R. Later Buchdahl published a paper about some Mathematical 
solutions with this form of action [2]. These theories could be thought of as a special kind 
of higher derivative gravitational theories. As in the case of the EinsteinHilbert action one 
could derive field equations in f (R) gravity in two different approaches, the so-called metric 
and Palatini approaches. But in f (R) actions, unlike the EinsteinHilbert action or its mod- 
ified version (one with a cosmological constant term), the field equations obtained by the 
two approaches are not the same in general. In what follows we will be interested only in 
metric f (R) theories of gravity in which connection is dependent on metric g^ u with respect 
to which, the action is varied. In ordinary GR there are not that many exact solutions of 
the field equations for a given symmetry. Being higher derivative theory it is not unexpected 
to find more exact solutions in f (R) gravity and this turns out to be the case for spheri- 
cally symmetric solutions [ 3] .C onstant curvature solutions in f(R) gravity was studied by 
Odintsov and collaborators [111] .Also there is a significant review of f(R) and other modified 
gravities 12] and viable models 3] in literatures . Former Looking for solutions with a 
different symmetry, we obtained static cylindrically symmetric vacuum solutions of the f 
(R) modified Einstein equations [1]. It is shown how one can reduce the set of equations 
into a single equation which could then be utilized to construct explicit solutions. Now we 
generalized field equations with the same symmetry For constant curvature solutions, using 
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the general form of a cylindrically symmetric solution in Weyl coordinates we find, among 
possible solutions, a generalized form of zero curvature spacetime as well as one new (non- 
zero curvature) solutions with one of their parameters chosen so that it is related naively 
to the cosmo logical constant. The stability of a non linear system of differential equations 
is a significance problem. The modern method for study the stability problem is theory of 
dynamical systems that first is used by Coley in cosmological models [l4j.Folowing Coley we 
applied this method for our model and obtained a constraint on the initial condition which 
under it , the system of non linear equations are stable. Our motivation in this paper is that 
we solved the non linear equations and obtained Tian solution and we prove that our exact 
solution is applicable to the exterior of a string {5]. All mathematical calculations is done in 
detail .Therefore this work complete our former paper and although contains some new fea- 
tures of constant curvature solutions in f(R) metric gravity in cylindrical symmetry. Unlike 
our previous article we do not limit ourselves to a specific gauge , rather we work in a 
general framework. Finally we present a prediction for the general form of this cylindrically 
symmetric solution in other modified f{G) ,f(R,G) models of gravity. 



II. FIELD EQUATIONS IN /(R) GRAVITY 

In this section we give a brief review of the field equations in f(R) gravity. The action 
for f(R) gravity is given by 

S = J(f(R) + £ m )V=gd A x (1) 

The field equation resulting from this action in the metric approach, i.e. assuming that the 
connection is the Levi-Civita connection and the variation is done with respect to the metric 
SV, is given by 

G^tu = R[iv — -Rg^u = Tj^, + jpTj^j (2) 

where the gravitational stress-energy tensor is 

T*, = * (^,(/(J2) - RF(R)) + F{R)^{g a » gf)v - g, u g afj )) (3) 

with F(R) = df(R)/dR and T™ the standard matter stress-energy tensor derived from 
the matter Lagrangian £ m in the action (1). The vacuum equations of motion, i.e. in the 
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absence of matter, are given by, 

F(R)R pu - lf{R)g^^V v F{R) + g^DF(R) = (4) 

Contraction of the above field equations gives the following relation between f(R) and its 
derivative F(R) 

F(R)R - 2f(R) + DF(R) = (5) 

which will be employed later both to simplify the field equations and to find the general 
form of the f(R) function. 

III. CYLINDRICALLY SYMMETRIC VACUUM SOLUTIONS 

Interested in the static cylindrically symmetric solutions of the vacuum field equations 
(4), we start with the general form of such a metric in the cylindrical coordinates (t, r, ip, z) 



given by 



Id; 



ds 2 = -A{r)dt 2 + dr 2 + B(r)d<p 2 + C{r)dz 2 (6) 



The corresponding scalar curvature is 



*= t(|4fj-^(f) 2 ) m 

(*7^J")=1 



in which ' = j- and Ai(r) = A(r), B(r),C(r) for i = 1...3 are the metric functions. Using 
equation (5), the modified Einstein equations become 

FR^u - V M V,F = \g^{FR - UF{R)) (8) 

As in the spherical case, since the metric only depends on the cylindrical radial coordinate 
r, one can view Eq. (8) as a set of differential equations for functions F(r), u(r), k(r) and 
w(r). In this case both sides are diagonal and hence we have four equations. Differentiating 
Eq. (5) with respect to r we have the extra consistency relation for F(r), 

RF-RF + 3(DF)' = (9) 

Any solution of Eq. (8) must satisfy this relation in order to be also a solution of the original 
modified Einstein's equations. From Eq. (8) it is found that 

FR m - V M V M F = 1 _ DF , R ,, (10) 



In other words the combination = - f - r ^~ v m v m- f (ynth fixed indices) is independent of the 
index \i and therefore = A v for all fi, v. This allows us to write the following independent 
field equations; 

1 A 1 R' C 1 R" /""' R' 

+ il A ( B + C ' - 2( 1T + C") + « B f + ( C ^ F = ° <"> 

+ 2 1^ A - 1T> + "B - < A » + C < A - 1? = ° < 12) 
1 r" 4// r?" nj a 1 u' a' nj 

+ 2 P( T " C » + « C » " "A + B ( A " » F = ° 

corresponding to A t = A r , A t = A^ and A 4 = A z respectively. Therefore, any set of 

functions F(r), A(r), B(r) and C(r) satisfying the above equations would be a solution 

of the modified Einstein field equations (8). Obviously these equations can not be solved 

without auxiliary conditions reducing the number of the unknown functions. In the following 

section we discuss the simple but important case of the solutions with constant curvature. 



(13) 



A. A brief discussion about cylindric constant curvature solution for couple real- 
istic F(R) models 

One reason for investigating the constant curvature solution in a typical modified mod- 
ified f(R) gravity is that , in order that the accelerating expansion in the present universe 
could be generated, we must consider that f(R) could be a small constant at present uni- 
verse [20(.In this class of models, the universe starts from the inflation driven by the effective 
cosmological constant at the early stage, where curvature is very large. As curvature be- 
comes smaller, the effective cosmological constant also becomes smaller. After that the 
radiation/matter dominates. When the density of the radiation and the matter becomes 
small and the curvature goes to the value Rq , there appears the small effective cosmological 
constant. Hence, the current cosmic expansion could start. In the epoch of the inflation, the 
curvature R = A could be large , we may identify /(oo) as the cosmological constant for the 
inflationary epoch and f(Ro) as that at the present accelerating era |2o| . 



B. Constant Curvature solutions in Cylindrical symmetry 

It is known that some of the vacuum constant curvature solutions in f(R) gravity are 
equivalent to vacuum solutions in Einstein theory with the same symmetry. For example 
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it is shown in [J] that in the spherically symmetric case the corresponding f(R) solutions 
include the Schwarzschild-de-Sitter space for a specific choice of one of the constants of 
integration. For cylindrical symmetry, in Einstein gravity, static vacuum solutions were 
found almost immediately after their spherical counter parts by Levi-Civita [s] but those 
with a cosmological constant have to wait another 60 years to be found by Linet 4( and 
independently by Tian [5[ in a non-Weyl coordinate system. Their solution reduces to that 
of a cosmic string in the limit r — > 7|, |8|. 

Looking for cylindrically symmetric solutions in f(R) gravity, here we consider the simple 
but physically important case of static constant curvature spacetimes. To do so, in the 
field equations (11), (12) and (13) taking R = constant, we arrive at the following set of 
equations: 

:<f) 2 +(^) 2 ) = o (M) 

+§<£-§)=<> («) 

These equations must be solved . One simple assumption is seeking power law solutions for 
metric functions i.e. A(r) oc r m% which leads to 3 simple algebraic equations for mj. which 
as we show in next section this solution and another solutions. In the rest of this section we 
show that this solution with another set of solutions correspond respectively to a zero and 
non-zero values of Ricci scalar R. 
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Case (1) : solution with R = 

From the equations (14-16) one could obviously arrange a solution of Ai(r) oc r mi in 
which the constants m« are given as follows; 

nil — m 3 — 0, m 2 = 2 (17) 
4 2 

mi = m 2 = -, m 3 = -- (18) 

It can be shown that these are Ricci flat solutions (i.e R = in (7)) .Indeed in the second 
case which mi = n%2 = |,J7i3 = — § the metric is written as: 

ds 2 = -A r A ^dt 2 + dr 2 + B r A / 3 dv 2 + C r 2 / 3 dz 2 (19) 
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by applying the following radial coordinate transformation 



4 
3' 



/ = -//• ; (20) 



and choosing the constants A = B = (|) 4//3 ,C = (|) 2 ^ 3 it transforms into the following 
historic metric 

ds 2 = -pdt 2 + /rV /2 (^ 2 + dp 2 ) + p 2 V] (21) 

it is Levi-Civita's static cylindrically symmetric solution which normally is written in the 
following general form with the constant m—\: 

ds 2 = -p 2m dt 2 + p~ 2m [p 2m \dp 2 + dz 2 ) + p 2 d^ 2 } (22) 



That it is obviously flat 



0. 



Case (2) : solution with R = constant / 

To find solutions of this type we introduce 3 new functions as A = e u , B = e v , C = e w . 
In terms of these functions the equations (11-13) become: 

u'(v' + w') - 2(v" + w") - (v' 2 + w' 2 ) = (23) 
2(u" - v") - (v' 2 - u' 2 ) + w'(u' -v')=0 (24) 
2(u" - w") + (u' 2 - w' 2 ) + v'(u' -w')=0 (25) 

Solving this equations is complicated. But in a special case u = v equation (24) is satisfied 
identically and 2 remaining equations(23,25) reduce to: 

u" — 2 U ' W ' ~^ w " ~^ 2 W ' 2 = ® 

// 1 / / // 1 /2 i /2 n /ri7\ 

U — —U W — W — -W +U =0 (27) 

By adding (26,27) we have: 

2u" + u' 2 - u'w' = (28) 
a straightforward integration gives: 

w = 2\nu' + u + c (29) 



in which c is a constant . After substituting this solution in (26) we have: 

/// 

u 

3u" + 2— = (30) 



this is a simple 2'nd order ODE w.r.t u'. Taking p = u' and rewriting the derivatives with 
respect to u we have: 

2 d P 1 2 /oi\ 

3 P chl = C -2 P (31) 

which has the following solution: 

p = ±^2c + Ae~ 3u / 2 (32) 
If we choose only positive sign we can write an explicit form for u, v, w: 

u=v= l in 4c cos2{ lS- r)] (33) 

3 2 A 3 

w = 2 ln[V2c. tan(- V2c.r)} + - ln[— cos 2 {-y/2c.r)] + c 2 (34) 

Without loss of generality we can set c 2 = or absorb it in the z coordinate. If we calculate 
Ricci scalar (7) with these metric functions we obtain R = 6c. So by comparing this with 
the solutions of the Einstein field equation in the presence of the cosmological constant i.e, 
R = 4A,we have: 

c = h (35) 
Finally we can write the general form of metric (6) as: 

ds 2 = -\j c cos 2 (^.r)] 2 /W + dr 2 + cos 2 (^.r)] 2 / 3 V + (36) 

[2A 2 csin 2 (3/4v / 2^r) cos^^v^r)]^ 2 

Now we introduce new parameter a = (^) 2 ^ 3 - Our metric solution converts to: 

ds 2 = acos 4 / 3 (^r)(V - dt 2 ) + dr 2 + (^Aa) 3 sin 2 (^r) cos 4 / 3 (^r)^ 2 (37) 

from the general form of Tian metric we have: 

ds 2 = -[tan/3(r + f)] 7l [sin2/3(r + f)] 2 / 3 rft 2 + [tan/3(r + r)] 72 [sin2/3(r + f)] 2 /W (38) 

+ [tan/3(r + f)] 73 [sin2/3(r + r)] 2/3 dy? 2 + dr 2 



Where f is an arbitrary real number, /3 = \/ ^ and 71, 72 73 are real numbers in the interval 

4 4] 
3' 3-1' 



1 ' 1 satisfying the algebraic equations 
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7i + 72 + 73 = 0, 7172 + 7273 + 7i73 = - g (39) 

Easily we observe that our solution is a special sub family of Tian solution with parameters 
71=73 = — |,7 2 = |,r = 0,a = . Also a change of coordinates (it is better to called 
them a congruity of non radial coordinates )is needed as: 

^(^W->(^f>* (40) 

Then we prove that in f(R) gravity the constant curvature solution in cylindrically sym- 
metric cases is only one member of the most generalized Tian family in GR, thus we have 
different kinds of cosmological constant in this scenario and not a fixed one as in GR. 



IV. MATHEMATICAL BEHAVIOR OF FIELD EQUATIONS AS AN AU- 
TONOMOUS DYNAMICAL SYSTEM 

In this section our goal is description stability of system of ODE (23-25) in modern 
language of dynamical systems [ijj]. A non linear system of differential equations X = 
fi(X),X = (X U X 2 , ....,X n ),X = % is called autonomous . A point P = {X%,X°, X%) 
is called a critical point if that be a solution for the system of algebraic equations: 

f i (X) = 0,i = l,2,...,n (41) 

any autonomous system can be spotted near critical point as a linear system. If we construct 
a matrix A with components Ay = ( 9 q^ )\p , and if the linear perturbation of system from 
critical point is denoted by \l/ = (^1, ^2, ^n), then the solution to linearized system 
^ = A^fr determines stability of the original non linear model. If as t — > 00 all perturbation 
functions — > then we say that the system is stabilized. For our purpose if in system 
(23-25) we introduce a new set of variables £ = u',rj = v',6 = w' and if we denote all of 
them in form i = 1, 2, 3 we can write this equations in form : 

= lititk - ^ k, 1,3, k = 1, 2, 3 (42) 
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it's related matrix A is: 

Aij = -tijk£,k — (43) 

Where in it e^j. is a Levi-Civita altering tensor and 8y is Kronecher delta function. The 
critical point of system is located at P = (£, ?7, #) = c(l, 1, 1) where in it c is an arbitrary 
constant. Thus it's value at critical point can be written explicitly by: 

A i=j = -c, A^j = X -c (44) 

the eigenvalues of matrix A are: 

Ai = 0,A 2i3 = -^ (45) 

Thus the general solution for system (42) in linearized aproximation is : 

# = £o + (?7o + 0Q)e-i r (46) 

For stability we observe that only if £ = the non linear system (42) in linear approximation 
is stable and consequently the non linear system begins stable. 



V. PHYSICAL RESULTS AND CONCLUSIONS 



In recent decades there has been more interests in Cosmic strings as an Cosmological 
predictions of String theory. We know that according to the gauge theories with spontaneous 
symmetry breaking the universe has a number of phase transitions since the big bang. In 
general, cosmic strings are topological structures which is produced meanwhile in these 
phase transitions. Later Grfinkle [15] has treated the string as a self-interacting massless 
minimally coupled scalar field coupled to a U(l) gauge field and has shown that there exists 
a class of static , cylindrically symmetric solutions to the field equations. Now following 
the Tian (5] we prove that our metric (37) is applicable to the exterior of a string. This 
is the main physics which is hidden behind mathematical equations and we elucidate it in 
the following paragraphes. Consider the exterior metric for a string. One class of exteriors 
has been discussed by Tian is a solution of equations in a vaccum or by imposing the boost 
symmetry. The boost symmetry means that in the general Tian metric we must have 

2 

7i = 72 = ±-, 
4 

73 = ±- 
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The result metric is: 

d 2 a~ 2/3 

ds 2 = [acos/3(r + f)] 4/3 (-dt 2 + dz 2 ) + — — — [sin 0(r + f)] 2 x [cos/3(r + r)r 2/3 dp 2 + dr 2 

p* 

If we apply the complex transformations t — > —itp, z it,(p —t z, this form accordance with 
the our metric (37) by choosing parameters f = , a = 3 j^ L , a = a 3//4 , d = /3a 1 / 3 . This proofs 
that our metric (37) is applicable to the exterior of a string. 

Again the motivation for considering material systems with cylindrical symmetry in asymp- 
totically AdS spacetimes comes also from the classical general relativity, it is well known 
that with A = , the asymptotes of cylindrically symmetric static spacetimes representing 
infinite sources is very different from that of spatially bounded static sources. The advantage 
of the studied spacetimes over the vaccum LC case is that for A < 0, far away from the 
symmetry axis, r = , they approach the anti de Sitter solution (AdS) unlike LC that does 
not approach the Minkowski spacetime. Lemos 16] found an interesting class of cylindrically 
symmetric 4D solutions with A < called black strings. These spacetimes describe the fields 
of charged, rotating strings that generally feature singularities and horizons and are asymp- 
totically anti de Sitter far away from the axis of symmetry. Static plane-symmetric perfect 
fluids occur also as subcases of the static cylindrically-symmetric solutions, There are sev- 
eral physically interesting classes of solutions admitting an Abelian group G(2) on non-null 
orbits; because of their importance and the large amount of relevant material, they will be di- 
vided into separate parts. The stationary axisymmetric fields have timelike group orbits T 2 . 
The classes with spacelike group orbits S 2 Contains time-dependent cylindrically-symmetric 
fields and their stationary subclasses, colliding plane waves with their typical dependence 
on the retarded advanced time, and inhomogeneous perfect fluid solutions. Known rigidly 
rotating axisymmetric perfect fluid solutions with more than two Killing vectors belong to 
the locally rotationally-symmetric space-times or to the homogeneous space-times or they 
are cylindrically-symmetric . The stationary cylindrically-symmetric dust solutions are con- 
tained in the general stationary axisymmetric dust solution (the Winicour solution) 
About generalization of our cylindrically symmetric result for other modified gravities es- 
pecially for f(G) and f(G, R) one could obtain new families of physically tolerable solu- 
tion. Former stationary cylindrically symmetric solutions to the five-dimensional Einstein 
and Einstein - Gauss - Bonnet equations has been investigated and two exact solutions 
which qualify as cosmic strings, one corresponding to an electrically charged cosmic string, 
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the other to an extended superconducting cosmic string surrounding a charged core were 
found[l8]. Thus naturally we expected that in any f(G) model of gravity by imposing the 
cylindrical symmetry one could be able to recover this rather old family of exact solutions. 
Note that in /(G)-gravity, there are no problems with the Newton law, instabilities and 
the anti-gravity regime |2l|.The analogue of constant curvature in /(G)-gravity is Go corre- 
sponds to the present value of the Gauss-Bonnet invariant. Technically when G = Gq and 
G = G+oo, f{G) becomes almost constant and can be regarded as the effective cosmological 
constant similar to f(R) models. The starting point for doing it is following from which has 
been done by Nojiri and Odintsov for spatially-flat FRW universe metrics in the context of 



modified Gauss-Bonnet f{G) and f(G,R) gravity[19j. As an landscape it seems t 



rat for a 
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2l| , the 



large number of choices of the function f(G), like that has been discussed in 
field equations has a non-trivial real solution (in which non of metric functions vanishes) and 
asymptotically converge to a static patch of deSitter universe. We remind it to the another 
work {22)]. 
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